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ABSTRACT 


A series of reacting flow tests of a solid fuel ramjet 
were conducted with different air inlet swirl. Cold flow 
measurements of the air flow at the exit plane of the tube- 
in-holé injector were used to determine the Swirl Number for 
each configuration. Regression rates of an HTPB type fuel at 
low air mass flux were found to increase with swirl vane 


angle (and Swirl Number) up to 30 degrees. 
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I. INTRODUCTION 


A. SOLID FUEL RAMJET 

Ramjet technology has taken three distinct approaches to 
meet the requirements of an ever-expanding mission responsi- 
bility placed upon weapon systems of today. These three 
concepts are the Liquid Fuel Ramjet (LFRJ), Solid Fuel Ram- 
jet (SFRJ) and the Ducted Rocket (DR). Each system has 
distinct advantages over the other. The LFRJ has the 
capability to operate over widely different levels of thrust 
and efficiency through use of variable fuel control. The 
SFRJ is lighter, more compact, and is much less complex 
than the LFRJ due to the absence of fuel tanks, fuel 
delivery and control mechanisms. The DR retains some of 
the benefits of the SFRJ, but must carry oxidizing materials 
within the fuel to maintain combustion. 

Because of the simplicity with which the SFRJ operates, 
and its relative compactness, it is desireable to develop 
a simple thrust control for the solid fuel ramjet motor. To 
understand how this problem may be approached, it is first 
necessary to have a working knowledge of the basics of ramjet 
technology. 

As stated earlier, the ramjet is a very simple engine. 
Since there is no compressor, the ramjet uses supersonic 


velocities to compress the air flow. Figure 1 depicts the 
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basic SFRJ. The three components of the basic SFRJ are the 
supersonic diffuser or inlet, the combustor, and the 
exhaust nozzle. 

As supersonic flows enter the SFRJ diffuser, several 
Significant events take place. Standing shocks are present 
at the air inlet. An increase in static temperature and 
decrease in flow velocity occurs through the shock waves. 
Combustion of the solid fuel increases the mass flow slightly 
and the temperature is again increased. This larger mass 
flow at higher temperature is then expanded through the 
exhaust nozzle. A detailed description is given by United 
Technologies Corporation, Chemical Systems Division [Ref. 1]. 
It is the increase (or decrease) of the fuel prolysis rate or 
regression rate (r) that determines the thrust and effects 
the efficiency of the SFRJ. 

There are many methods that can be used to affect the 
regression rate of a solid fuel. One is to change the fuel 
composition. Another is to change the environmental effects 
on the fuel, such as’ the chamber pressure (Pc), inlet air 
temperature, and the mass flux through the port of the fuel 
grain. A typical regression rate formula is given by Boaz 


and Netzer [Ref. 2] 


for motors not uSing bypass air where 
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c = constant 
Pe = chamber pressure 
G = mass flow per unit area in the fuel port 


Ti = inlet air temperature 


The SFRJ is self-throttling through the effect of G on 
ae However, the self-throttling is not sufficient for many 
missions. Alternate means for varying r are therefore 
desired. Changing G can be accomplished by either changing 
m (such as in bypass designs) or by changing the flow area. 
Variable bypass is one possible means for controlling r. 

Ko [Ref. 3] investigated the effects of inlet air swirl 
on the regression rate of the solid fuel grain. Swirl 
imparts an angular component of mass flux into the flow field 
which can effect the residence time of the flow by increasing 
the effective length of the combusion section of the ramjet 
motor. It can also affect the convective heat transfer to 
the fuel surface and the mixing rates of fuel and oxygen. It 
was found that swirl may be an effective means for regression 


rate control, although only limited data were obtained. 


B. SWIRL NUMBER 

The introduction of swirl into a flow causes difficulty 
in characterization of the flow field. There is no longer 
only the axial component of the flow but now radial and 
tangential components as well. The Swirl Number (S) is one 


way to characterize the flow. Lilley and others [Ref. 4] 


ne 


have done considerable work with swirling flows and have set 
forth the non-dimensionalized swirl number (for specific 


assumptions about the inlet swirl behavior) as 





G 
an Oa (2) 
G. (d/2) 
where 2S 
G. = tangential flux of swirl momentum 
Gx = axial flux of momentum 
a/2 = equivalent inlet radius 


This thesis investigated the effects of varying the inlet 
air flow angle on the combustion process of a non-bypass 
SFRJ. This was accomplished by using tube-in-hole injectors 
with vanes of 0, 15, 30, and 45 degrees. Cold flow 
measurements were made to determine inlet flow angles, and 
reacting flow tests were conducted to determine the effects 


on fuel regression rate. 
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II. DESCRIPTION OF APPARATUS 


A. RAMJET MOTOR 

The ramjet motor used in this investigation was identical 
with one used in earlier invesigations at he Naval 
Postgratduate School [Refs. 2,3]. Figure 2 is a schematic of 
the head-end assembly, step insert section, fuel grain and 
nozZle. 

The head-end assembly used a wedge to turn the primary 
inlet air from 90 degree plenum dumps. The igniter torch and 
and ignition gas were injected into the recirculation zone 
just aft of the step insert. 

The step insert section acted as the air inlet, which 
was modified with a tube-in-hole injector. The step height 
allowed for sudden expansion of the flow into the combustion 
section as depicted in Figure 3. This sudden expansion 
allows a flame to be stabilize behind the step by creating a 
recirculation zone. The swirl elements and tube were 
fabricated from stainless steel. Each swirl blade was a 
cambered airfoil attached at the hub with silver solder. The 
vane tips were ground for minimum clearance at the tube inner 
wall, and then silver soldered along the downstream side of 
the blades. Figure 4 shows the tubes with the swirlers used 


in this experiment. Tubes with 15, 30, and 45 deg. (from 
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Figure 4. 15, 30, and 45 Deqree Swirl 
Blades Mounted Inside Tubes 
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the axial centerline) swirl vane elements and a tube without 
a swirler were used. 

CSD-18818 (an HTPB based system) was the solid fuel used 
for the hot runs. The average length and internal diameter 
of the fuel grains used were 12.0 and 1.7 inches 
respectively. The fuel grain was placed between the head 
assembly and an aft mixing chamber. The latter contained the 
exhaust nozzle and was connected with steel rods which bolted 
to the head-end assembly. The entire motor was mounted on 
the thrust stand. Figure 5 is a photograph of the thrust 
stand a the test cell. The control panel was used to 
initiate the runs and the Hewlett-Packard 9836 Computer and 
3054A Data Acquisition/Control system sequenced events during 


the runs. 
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Figure 5. Experimental Thrust Stand and SFRJ 
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III. EXPERIMENTAL PROCEDURES 


A. CALIBRATION 

Calibration of pressure transducers was conducted 
at the beginning of each test day utilizing a dead-weight 
tester. Table 1 is a list of the transducers calibrated and 


the calibration pressures/load. 


TABLE 1 


TRANSDUCER CALIBRATION PRESSURES/LOAD 


Pa 500 psi 
Pe 100 psi 
Ph 100 psi 
Phf 400 psi 
Pho 500 psi 
Pif 250 psi 
F 100 l1bf 


B. DATA EXTRACTION 

A Honeywell 1508B Visicorder was used to obtain analog 
readings from pressure and thrust transducers. Chromel- 
alumel thermocouples were used for measurements of all gas 
temperatures. A Hewlett-Packard 9836 Computer and 3054A 
Automatic Data Acquisition/Control Stytem were utilized to 


extract digital information. 
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C. HOT FIRING 

Fuel grain length (L), average port diameter (Dp) , and 
initial weight (W) were measured for each run. The inlet was 
then mounted on the ramjet motor and preparations were made 
for the ensuing runs. 

Air and ignition gas pressures were calculated using the 


choked ‘flow formula (assuming Cd = 0.97) 


(3) 





(3a: 
Ca PL 4 


a rm 
(R Ty) : 
Pressures were set using dome loaders and mass flow rates 
were set using the computer system as the controlling device. 
The computer was then readied for data acquisiton and each 
test was initiated from the control panel. Table 2 lists the 
inlet geometry and test condintions for each run. 

Upon completion of each run, the fuel grain was weighed 
and weight loss recorded. Actual motor burn time (tb), air 
time (ta), ignition time (ti), pre-ignition chamber pres- 
sure and thrust (Pa and Fa), and average chamber pressure 
and thrust (Pc and F) were obtained from the traces recorded 
on the Visicorder. The change in fuel grain weight is be 
expressed as 

AW = 7 (dz - d?)Lo, (4) 


Solving for the average final port diameter, 


A ee Ne (Sy 
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TABLE 2 


INLET GEOMETRY AND TEST CONDITIONS 


Run Tube Length Swirl Angle ta,ti,tb,tp 
(in. ) (deg. ) (sec.)** 
1 25 e) 6,1, 6,4 
2 Zo 0 6,1, 6,4 
3 PAP: 15 6,1, 6,4 
4°% 245 15 6,1, 6,4 
5 25 30 6,1, 6,4 
6 25 30 6,1, 6,4 
7 25 45 6,1, 6,4 
8 255 45 6,1, 6,4 

Remarks: 


Unsuccessful first two attempts. Cooled motor, 
reweighed grain, cleaned ignition gas choke and 
torch. Third try successful. 


pre-ignition airflow time, ignition time, burn time and 
purge time respectively 
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and substituting into 


(a- - d.) 
a eee (6) 


2% 


yeilds the average fuel regression rate. 


ry I 
| 


Average values for m(air), m(O,), m(C H,). T. and 
m(ign fuel) were then calculated from the digital printout 


for the run. The fuel mass flow rate was calculated by 


aatle 


— 4 
The Naval Weapons Center (NWC) China Lake, Ca., 

Propellant Evaluation Program (PEPCODE) was used to obtain 
theoretical values for the adiabatic combustion temperature 
and combustion gas properties, Y and R, using Pc, Ti, 
heats of formation, and elemental composition by weight of 
the fuel grain and injected gases used during each run. Ko 
(Ref. 3] gives a complete explanation of the procedures used 


to obtain the values needed for data reduction. 


D. COLD FLOW PROCEDURES 

A hemispherical, three-hole stagnation pressure probe 
was used to measure flow velocities at the exit plane of the 
tube-in-hole injector. The probe was connected by tubing to 
a pressure transducer. The transducer in turn, was connected 
through a controller to a voltmeter where stagnation pres- 
sure readings in millivolts were displayed. A Scanivalve 


allowed for cycling through the three pressure readings to 
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find the angle at which the flow was eminating from the 
tube-in-hole injector exit plane. This was accomplished by 
finding the maximum value for the center port and nulling the 
left and right values. The velocities between the tube and 
the inlet wall were also measured. 

The probe could be rotated about its own axis. It was 
mounted on a vernier translation stand to allow for movement 
across the flow in 0.05 in. increments. The probe itself 
could be translated on its own axis as well as rotated. 
This facilitated taking measurements in the y-z plane with 
rotation about the y-axis. Figure 7 is a photograph of the 
probe mounted on the thrust stand in this configuration. 
With the angle of flow taken at each grid point, the probe 
was then rotated from the y-axis to the z-axis. Figure 
8 depicts the probe head and the degrees of freedom in which 
it could be moved. 

The pressure readings in millivolts and the angles at 
which the the maximum values occurred were then input into a 
program that caluclated x, y, and z components of the flow 
velocity, and transformed them from the Cartesian coordinate 
system to the cylindrical coordinate system. This produced 
axial and tangential velocity components (assuming that the 
net radial component was equal to zero). The computer 
program and examples of values calculated are given in 


Appendix A. 
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Figure 7. Photograph of Pressure Probe on Stand 
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Figure 8. Pressure Probe Head and Degrees of Freedom 


Upon completion of characterization of the flow at each 
point in the grid, the amount of swirl in the flow was then 
characterized using Swirl Number (S) as previously stated 


in equation 2. The tangential flux of momentum is 


3 
G, = +? UW (4/2) (8) 


and thé axial flux of momentum.is 


_™ pu_ (a/2)* (1 - (g/2)*) (9) 
2 


G oO 


x 
Making the same assumptions as Lilley [Ref. 4] of solid body 
rotation (a flat upstream axial velocity profile, and an 


increasing Swirl velocity from 0 at the hub to w,, at the 


outer edge), the analysis yeilds the Swirl Number given as 
= eee) 2 ae 
2 = T= (g/2)°) rer 
where 
“mo 
dS (11) 
Mo 


is the ratio of the maximum tangential velocity to axial 
velocity at the exit plane of the tube-in-hole injector. 

In an attempt to more accurately characterize the swirl 
in the present apparatus, a maSsS averaged momentum ratio was 


calculated by taking the mass average tangential velocity 


_tmwo (12) 


= 
lI 


Ze 


and the mass average axial velocity is given by 


= cmu 
Sr ii (13) 


Dividing equation 12 by equation 13 yeilds S(avg). 
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IV. RESULTS AND DISCUSSION 


A. INITIAL TESTS 

A series of initial hot runs was completed to test the 
integrity of the tube-in-hole injector design. During this 
phase it was found that initially designed tubes, which 
extended 0.25 inch beyond the aft end of the inlet, did not 
allow ignition of the motor. It was suspected that swirling 
flow was interferring with the recirculation zone, destroy- 
Ong the flame stabilization area. The tubes were lengthened 
so that the swirl dump was 1.5 inches beyond the inlet step, 
which solved the problem. This length was considerably less 
than the reattachment length of approximately 3.0 inches. 

Burn time for the motor was initially set for 10 seconds, 
which almost resulted in for burn-through to the wall. Burn 
time was therefore reduced to 8 seconds, which worked well 
for non-swirl test runs. Burn-through occurred at 8 secs. 
with 30 degrees of swirl, so the burn time was further re- 


duced to 6 seconds and the remaining hot runs were conducted. 


B. COLD FLOW ANALYSIS 

A summary of results of the cold flow studies is given in 
Table 3. There was an increase in the difference between the 
maximum and mass averaged axial velocities from 8% for no 


Swirl, to 43% for 45 degrees of swirl. The same was true 
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for the difference between the maximum and mass averaged 
tangential velocities, with 21% difference at 15 degrees and 
43% difference at 45 degrees. This could have heen a 
function of the increased blockage within the tube. The 
velocities between the tube and inlet wall increased as the 
vane angle increased, which was an indication that blockage 
was occurring. The Swirl number, as calculated by Equation 
10 (based on solid body rotation), increased with vane angle 
to 30 degrees and then decreased. The mass averaged value 
(Equation 12 divided by Equation 13) increased continually 


with vane angle. 


TABLE 3 


COLD FLOW RESULTS 


Blade angle 0 iS 30 45 
(deg) 

u(max ) 475 616 759 820 
(ft/sec) 

w(max) 0 309 S21 525 
(ft/sec) 

u(avg) 434 509 560 462 
(ft/sec) 

u(outer flow) 384 393 503 5a9 
(ft/sec) 

w(avg) 0 243 345 296 
(ft/sec) 

S(avg) One 0.478 0.616 0.640 
S On0 0.335 0.649 0.544 


a2 


C. FUEL GRAIN PHYSICAL CHARACTERISTICS 

The fuel grains used without swirl exhibited smooth 
burning surfaces with the port diameter slowly increasing to 
a point approximately 2.0 inches beyond the step. There- 
after, the port diameter was approximately constant. The 
grains used for the 15 degree swirl showed the same charact- 
eristi@s. However, four areas that corresponded with the 
four vanes of the swirl element exhibited slightly indented 
burn regions extending from 2 inches beyond the step to 4 
inches, slowly blending into a concentric surface. The 30 
degree grains showed a marked increase in the indented burn 
regions, and they extended to 5 inches heyond the step. The 
45 degree grains were even more pronounced in the indented 
burned region but not as elongated. The angles of the 
indented burned regions corresponded to the respective swirl] 
vane angles. Table 4 lists the physical characteristics of 


fuel grains pertaining to the tests. 


D. COMBUSTION PROPERTIES 

Measured quantities for each hot run are listed in Table 
5. A plot of the regression rate vs swirl vane angle is 
given in Figure 9 and vs’ S and S(avg) in Figure 10. | 
Regression rate increased rapidly for small amounts (15 
degree vane angle) of swirl. Further increases in swirl vane 
angles (or Swirl Number) had much less effect on r. The 


Swirl Number (Equation 10) correllated reasonably well with 
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regression rate. 


45 degree vane, 


A typical value of G (mass flux) in an SFRJ is 0.6 


lbm/in’-sec without bypass. 


yeild a fuel port G of 0.3. 


25% reduction in r could be achieved. 


the present configuration, 


this change in r. 


However, 


With rv Ge* 


A 30 degree swirl, 


it does appear that variable 


S(avg) appeared to be in question for the 


perhaps due to the large amount of blockage. 


Fifty percent bypass would then 


an approximately 


in 


would yeild only about half of 


swirl can be an effectiv method of regression rate control in 


the SFRJ. 
Run # Vane 
Angle 
(deg) 
1 0 
2 0 
3 3, 
4 L5 
5 30 
6 30 
7 45 
8 45 


L 


(in) 


122010 
12.000 
12.008 
12.014 
12,050 
L2 2060 
112929 


Lie 91D 


TABLE 4 


Dp 


(in) 


NE 310. 
1.749 
1.747 
Le oe 
li. 24 
1.740 
ET S18) 


1.746 


On 22 


0792 
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FUEL GRAIN PHYSICAL CHARACTERISTICS 


Dth 


(in) 


0.943 
0.943 
0.943 
0.943 
0.943 
0.943 
0.943 


0.943 


(sq. 


Ath 


0.6483 
0.6604 
0.6554 
0.67.35 
05,6597 
0.6525 
025721 


0.6648 


in) 


AW 


(gm) 


137 
200 
219 
PA ss 
2335 
230 
22 


220 


ti 
(sec) 


OGaZ 
0.84 
OL 86 
O.7S 
O273 
0.80 
0.79 


On87 


m 
air 


(lbm/sec) 


Oe 503 
O25.06 
0.495 
0-510 
0.508 
0.53 
oO; 303 


0.487 


tb 
(sec) 


| 


m 
fuel 
020598 
0.0608 
Om0G7 2 
0.0680 
0: 0713 
0.0702 
0.0602 


0.0667 


TABLE 5 


MEASURED COMBUSTION PROPERTIES 


ta Pa r EC 
(sec) (psia) (in/sec) (psia) 
6.0 5023 0.0236 POOR 2 
Gat 49.9 0.0240 O76 
6.0 49.5 0.0267 HOS. 7 
6.1 50.3 0.0266 107.8 
oon C) a eral 0.0277 108.6 
6.1 on Meg | 0.0274 HOO 3 
6.8 SOG 0.0263 HO728 
6.0 49.1 0.0276 HOB 
m m mn 
25 S35 if 
OPS LOLS, 0.0046 0.0033 
0} (Sh 0: 0.0047 O20535 
OOo} 0.0047 0.0034 
0.0193 0.0048 0.0034 
O,Oro1 0.0047 0.0034 
0.0194 0.0049 OZ e033 
0.0194 0.0048 0.0034 
O20192 0.0049 0.0036 
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V. CONCLUSIONS AND RECOMMENDATIONS 


It has been found that the regression rate increases 
Significantly with small amounts of swirl in the flow of 
inlet air, but that larger amounts of swirl yield less 
effect... The Swirl Number S (as given by Equation 10), at 
the low mass flux values of approximately 0.25 which were 
used in this investigation, correlated reasonably well with 
regression rate. 

It is recommended that the effects of Swirl Number on 
regression rate be investigated with increasing values of 
mass flux to 1.0 lbm/in -sec and with metalized fuels. 
Improvement of the 45 degree swirl element design should be 
also attempted, as well as an investigation of the proper- 
ties of combustion with swirl angles between 30 and 45 


degrees and beyond. 
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APPENOIX A 
PROGRAM FOR CCLOD FLOW ANALYSIS OF SWIRLING INLET AIRFLOWS 


SWIRL NUMBER PROGRAM 
LT.» WILLIAM He. CAMPBELL, JR. 
NAVAL POSTGRADUATE SCHOOL, MONTEREY 


¥ 


THIS PROG RAM IS USED TO COMPUTE THE AVERAGE SWIRL NUMBER OF A 
SWIKLING FLOW OF A TUBE IN HOLE INJECTOR FCR A SOLID RCCKET 
RAMJET AND TRANSFORM THE CARTESIAN VELOCITY COMPCNENTS INTO 
CYLINDRICAL VELOCITY COMPONENTS. 


DeFINITIONS 
TIT) = INLET AIR TEMPERATURE 
SONVEL (1) = SONIC VELOCITY 
DENS(1I) = AIR DENSITY 
YC 1K) = Y-CARTESTAN COORDINATE OF THE EXIT PLANE 
Z( Jel) = Z-CARTESTAN CCORDINATE OF THE EXIT PLANE 
RL(JsK) = RADIUS TO POINT YIN),Z(M) 
RY (Is JK) = PRESSURE VECTOR IN MILLIVOLTS ABOUT THE Y-AXIS 
RZCIy JK) = PRESSURE VECTOR IN MILLIVOLTS ABOUT THE Z-AXIS 


ALFACI »J5K) 
BETACT »J9K) 
Al( I,J 2K) 
A2(TsJ 9K) 


ANGLE MEASURED OF RY(NoM), POSITIVE ABOVE THE X=) 
ANGLE MEASURED OF RZ(NgM), POSITIVE LEFT QF Y=) 
X COMPONENT JF RYN) 
X COMPONENT DF RZONgM) 


5. 


B( I,JeK) 
C(IyJ-K) 
XCOMP(19J3sK)= 
RTGT(I Jo) 


RXZ(IleJSeK) 


a 


KXY(L a Sek) 


GAMMA(1sJ3K)= 
VELTA( I+J,K)= 
PHI(I eJsK) = 
TCOMP(IyJsK)= 


RCOMP( IsJ 7K) 
PTX(LaJeK) 
PTT (lL »JeK) 
UCI sJ9K) 

W( TJ, K) 
SUM1 = 


SUM2 = 
Sis = 


Y COMPONENT OF RZ(NeM) 

Z COMPONENT OF RY(N,%) 

TCTAL X COMPONENT OF PRESSURE IN MILLIVOLTS 
TOTAL VECTOR OF PRESSURE IN MILLIVOLTS 

IN MILLIVOLTS 

TOTAL VECTOR COMPONENT OF PRESSURE IN THE X-Z 
PLANE IN AILLELTIVOLTS 

TOTAL VECTGR CCMPONENT GF PRESSURE IN THE X-Y 


PLANE IN MILLIVOLTS 


ANGLE BETWEEN Y=0 AND RYZ 

ANGLE BETWEEN Y=0 AND Rl 

ANGLE S3ETWEEN X=9 AND RXY_, SWIRL ANGLE UF FLOW 
TANGENTIAL COMPONENT GF STAGNATION PRESSURE 

IN MILLIVOLTS 

RADIAL COMPONENT OF STAGNATION PRESSURE IN MILLIVULTS 
X-COMPONENT OF STAGNATION PRESSURE IN PSI 
TANGENTIAL COMPONENT OF STAGNATION PRESSURE IN PSI 
AXIAL VELOCITY COMPONENT 

TANGENTIAL VELGCITY COMPONENT 

SUMATION GF THE GRID REPRESENTATIVE AREA 
MULTIPLIED BY THE GRID POINT TANGENTIAL 

VELOCITY W(I,J+K) 

SUMATION OF THE GRID REPRESENTATIVE AREAS 

SUMATION OF THE GRID REPRESENTATIVE AREA MULT. 

BY THE oRIC POINT AXIAL VELOCITY JU(I,JsK) 
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UMAX(T) = THE MAXIMUM AXIAL VELOCITY FOR THAT RESPECTIVE 
SWIRL VANE ANGLE 

WMAXCT) = THE MAXIMUM TANGENTIAL VELOCITY FOR THAT 
RESPECTIVE SWIRL VANE ANGLE 

AW(T) = THE AVERAGE TANGENTIAL VELOCITY FGR THAT 
RESPECTIVE SWIRL VANE ANGLE 

AU(T) = THE AVERAGE AXIAL VELOCITY FOR THAT RESPECTIVE 

fs SWIRL VANE ANGLE 

SAV(T) = THE RATIO OF THE AW(I) TO AU(TI) 

SMAX(I ) = THE SWIRL NUMBER ACCIRDING TO LILLEY ET AL 

Gel = THE RATIO OF THE WMAX(I) TO UMAX(1) 


DELCARATION OF MATRICES FOR RESPECTIVE COMPONENTS CF THE PROGRAM 


REAL Y(1s8)5 Z(8eld, TI(4)-9 THETA(4), SONVEL(4), DENS (4),AW(4) , 


RY (49898) » ALFA(4 8 e8)_ RZ149898), BETA(4s848),AU(4), 
XCOMP (4878), R1(898)y RYZ(49898)_y RTOT (49898) eSUM1(4), 
PH1(4 9898) 4 RCOMP(498s8)5 PTX(4,893) 9 UlG%s8e8)_, SUM2(4) 
PTT(45 898) 9 W(14e898)4 G14), AREA(8,3) 4 SUM3(4)5 SUM4(4), 
Al( 41858) 5 A2(4%,898)9 Bl4rde8)y C141 898) eRXZ(448 98) 
RXV¥(498e8) 9 DELTA(4,898) 4 GAMMA(4 9858) 9 TCOMP(4,8,8), 
WMAX(4)—9 UMAX(4) ¢ SAV(4)_ SAAX(4) 


INCREMENTATION OF THE Y AND Z AXIS 


DO Ll N=1,8 
READ(5,10) Yl eN) 
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L(Ny 1) =Y(L_N) 
CONT INUE 
00 81 J=1y8 
READ(5 +400) (AREA(JsK)y K=1,98) 
CONTINUE - 
COMPUTATION OF SONIC VELOCITY,DENSITY, ANO INPUT OF VANE ANGLE. 
DO 2 1=194 
READ(5,20) TI(I), THETA(I) 
SONVEL (1)=SQRT(1 04%32.2*53.3¥*TI(1)) 
DENS(I ) =14-7/ (53 e34TI(1)) 
CONT INUE 


THIS PORTION OF THe PROGRAM INPUTS AND CUTPUTS DATA RECORDED 
OURING THE COLD FLOW ANALYSIS OF DIFFERENT SwIRLER VANE 
GEOMETRIES-TH VALUE (8+8) FOR EACH ANGLE IS THE COMPONENT 
MEASURED BETWEEN THE INLET WALL AND THE TUBE-IN-HOLE INJECTCR 
OUTER WALL. 


DO 4 [=1+4 
SUM1 (1) =0 
SUM2 (1) =0 
SUM3 (1) =9 
SUM4 (1) =0 
UMAX (1) =0 
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WMAX (1) =0 
DO 5 K=1,8 
D0 6 J=1,8 
READ( 5,30) RY(IsJsK)s ALFACTyJeK)5 RZGTeSeK)y BETACT yJyK) 

ALC IT sJeKI=ERYC 19 Je K)¥*COS (ALFACI yJsK) ¥32141/180) 
A2(TeJesKI=KZ¢0 19 Js K) *COS (BETACI » Je K) *32141/180) 
B(IyJeK) =RZC0 [yp JeK)¥SIN(BETACT »JpK)*32141/180) 
CCl ep JeKI=RV (1 pS eK) *¥SIN(ALFACT, JyK) *32141/180) 
XCOMP(I yp Je KI=AL( Ty Se KItAZ LI 9d9 K) 
R1L( JSeK) =SQRT( VOLK) **2+205,1) *¥2) 
RYZ (Lede KI =SORTI(C(C1L, JK) **¥24+B( Td 9K )*¥2) 
RXZ(IpJSeK) =SQRTIC( 1, Jp K) ¥*®24+XCOMPU I Sy K)**2) 
RXY (I gd 9K) =SQRT (XCOMP (1 oJ K)**¥24+B (TI 9 dy K) *¥*2) 
RTOT (Ts J9K) =SQRT(XCOMP (Is JeK)* *24+B( 1,5, K)**24+C (1 ,5,K) ¥¥2) 
IFC CY(L,K)-EQ-0) eAND.(Z0J91)2-EQ.0)) THEN 
DELTA(IsJK) =0 
ELSE IF ((Z0Js1)-GE.9.05) «AND. (Y¥(15K).6Q.0)) THEN 
DELTA(L,JS»K) =90 
EE SE 
DEL TACT» JSJ 9K) =(ATAN(Z( Ses LI/SYOL »K))) ¥180/3.14 
END IF 
IFC BCI sJsK)-EQ20) THEN 
GAMMA(I,JsK) =0 
ES E 
GAMMACT sJ eK) =CATANICOL, J eKI/BC 1, 5,K)))*1890/3.14 
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END IF 
IF (XCOMP(IyJdeK)eEQ.0) THEN 
PHI(1eJ9K)=0 
ELSE 
PHI(LeJeK)=C ATAN(B( 19d 9K) /XCOMP( Ty JK) 1)*180/3014 
END IF 
TCOMP( Ly JoK )=RYZ (19 Je K) #C9S ( (99 = (DELTA Lyd 9K) +GAMMA(T oJ 9K) )) 
+*3.141/180) | 
RCOMP(IpJeK)=RYZ( 19 de K) ¥SIN((90-(DEL TAC 19d 9K) +GAMMACT 9deK) )) 
+*3. 141/180) 
PTX ( LeJeK) =XCOMP( 1 9d 2K 720155 
U (Ls JS eK) =SQRT( 2402 6 7#TI ALPE (5#( C(PTX (19d 9K) 14. 7) +1) #628571 
+#1))) 
PTT (IeJdeK) =TCCMP( I pd eK 1/0155 
WIL pS oK)=SQRT (24020 THTIC IT) RI SE CCC PTT (Led oK IS 1407) FL R* .28571 
+-1))) 
SUM1 (1) =SUMLO I) #(AREA(J 9K) #WUT Je K)) 
SUM2(1)=SUM2( 1) #AREA (JK) 
SUM 3 (1) =SUM3(1)+( AREA(S, K) #UCL 9d eK) D 
IF (U(L9deK).GTeUMAX(1)) THEN 
UMAX( 1) =U (I) dy K) 
ELSE 
UMAX( 1) =UMAX(T) 
END IF 
IF (WOL9deK)eGTeWMAX (1)) THEN 
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WMAX( I) =W( Ty Jo K) 
ELSE 
WMAX(I) =WMAX (1) 
END IF 
CONTINUE 
CONTINUE 
AW(1)=SUML(I)/SUM2(1) 
AU( 1)="SUM3 (1) /SUM2 (1) 
SAV (I) =AW(I)/AUCT) 
G(1)=wMAX(1)/UMAX(1) 
SMAX(L) =(G01) 72) / 01-6 (1) #*2) 
CONTINJE 
WRITE(64410) (UMAX( I), L=194) 
WRITE(6 9411) (WMAX(I), [=194) 
WRITE(6 5412) (AW(I), 1=194) 
WRITE(6 4413) (AUC T), [=194) 
WRITE(6,414) (SAV(I), I[=1)4) 
WRITE(65415) (SMAX( I), I=1 94) 
DO 8 I=1,4 
WRITE( 6,40) THETA(TI) 
WRITE( 6,50) (Y¥CleL), L=ly 8) 
WRITE( 6,55) 
DO 9 J=1,8 
WRITE (6,60) ZlJol)s (UCT JeK) K=193) 
CONTINUE 
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Ll 
l2 


ee 
14 


15 
16 


CONT INJE 
DO 12 T=l+4 
WRITE@6570) THETA(T) 


WRITE(6,50) CY(1L»sL)»L=1,8) 


WRI TE( 6555) 
00 11 J=1,8 


WRITE (6960) ZlJolds (Wil sd sK),K=1,8) 


CONTINUE 
CONTINUE 
DO 14 1=194 
WRITE( 6,80) THETA(I) 


ARITE (6,50) (Y¥( 1,L) ,L=1,8) 


WRITE (6,55) 
DO 13 J=1,8 
WRITE(6e90) ZiJel)ds 
CONT TINUE 
CONTINUE 
DO 16 1[=1ls% 
WRITE( 6,80) THETA(T) 


(RY (I eJeK), K=1,8) 


WRITE(6550) (Y¥(1eL),L=1,8) 


wRITE( 6,55) 

VQ 15 J=198 

WRITE (6,120) ZU Jol), 
CONTINUE 
CONTINUE 


(PHI (I sJe2K) 9 K=1,8) 
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DO 18 1=1e4 
WRITE (69100) THETA(I1) 
WRITE (6,50) (Y(€15L),L=1,8) 
WRITE (6,55) 
DO 17 J=178 
WRITE(6990) Zl Jol)_9 (RZGLeJSeK)y K=1,y 8) 
17 CONTINUE 
18 CONT INUE 
DO 22 1 = 1le4 
WRITE (6,110) THETA(1) 
WRITE(6950) (Y(1yL) sL=198) 
WRITE (6955) 
DO 21 J=1;8 
WRITE (6,120) Z{Jdol)de (ALFACI+JSeK)_ K=1,8) 
21 CONT I NUE 
22 CONTINUE 
DO 24 I=l+4 
WRITE (65130) THETA(I) 
WRITE(6,50) (YC lsl), L=1,8) 
WRITE (6,55) 
DO 23 J=1,8 
WRITE(6,120) Zldol)de (B3ETACT oJ K)e K=128) 
23 CONT INUE 
24 CONT INUE 
ema 2 [=1,4 
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WRITE (6,140) THETA(I) 
WRITE (6950) (Y(1l9L), L=1,8) 
WRITE (6,55) 
DO 31 J=1+8 
WRITE(69150) ZlJdol)y (Cll gd9K)y K=1,8) 
31 CONTINUE 
32 CONTINUE 
DO 34 L=194 
WRITE (6,160) THETA(I) 
WRITE (6950) (YC lel), L=1,8) 
WRITE (6955) 
D0 33 J=198 
WRITE(69150) ZlJdol) 9 (XCOMP(IyJyK)y K=1,8) 
33 CONTINUE 
34 CONT INUE 
DO 36 I=1ls4 
WRITE (69170) THETA(L) 
WRITE (6550) (Y(19L) pL=158) 
WRITE (6,55) 
DO 35 J=1,8 
WRITE(69150) Z(Jol)y (RTOT(IL9JdyK)y K=1_8) 
35 CONTINUE 
36 CONTINUE 
DD 38 1=194 
WRITE (69180) THETA(1) 
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Sie | 
38 


41 
42 


43 
44 


WRITE (6950) (YC1,L), L=1,3) 
WRITE (6955) 
D0 37 J=1.8 
MRITEC6y 120) 2tdyl), (PHI¢l,JeK), K=1y8) 
CONT INUE 
CONTINUE 
00 42 1=1)4 
WRITE (6,190) THETA(I) 
MITE (6,50) CY( lel), L=1.8) 
WRITE (655) 
DO 41 J=118 
MRITE (6,150) Zde1)> (Al€l+JeK) y K=1098) 
CONTINUE 
CONT INUE 
D0 44 T=1s4 
WRITE (6,200) THETA(T) 
Mette (6450) (Y(l,L), L=1,8) 
WRIT E(6,55) 
DO 43 J=1,8 
MHEG, 150) ZtJells (AZ(IedeK) > K=1.8) 
CONTINUE 
CONT INJE 
DO 46 I=ly4 
WRITE(6,210) THETA(T) 
WRITE( 6950) (Y(LyLd, L=1,8) 


49 


45 
46 


47 
48 


51 
De 


WRITE (6955) 
DO 45 J=1,8 
WRITE (69150) ZlJol), (Bly JdeK)y K=1,8) 
CONTINUE 
CONTINUE 
00 48 1=1.4 
WRITE (6,220) THETA(T) 
ARITE (6950) (YO Let)» t=l 8) 
WRITE (6955) 
DO 47 J=1s8 
WRITE(69150) Zldol)_9 (RYZ(TedeK)_ K=1,8) 
CONTINUE 
CONTINUE 
00 52 [=ls4 
WRITE (69230) THETA(I) 
WRITEC 6,50) (YC let), L=L,8) 
WRITE (6,55) 
D0 51L J=1,8 
WRITE (69150) Z(del)s (RXZ41 9 deK) » K=1_8) 
CONTINUE 
CONTINUE 
039 54 [=ly4 
WRITE (64240) THETA(I) 
ARITE(6950) (VC lel), L=L98) 
WRITE (6955) 
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00 53 J=158 
WRITE (6,150) Zldel)» (RXYCI9J9K)y K=1¢8) 
53 CONTINUE 
54 CONT INUE 
D0 57 [=l»4 
WRITE (6,250) THETA(T) 
WRITE (6950) (YC lesb)» L=1,8) 
WRITE (6555) 
DO 56 J=118 
MRIDE(G,120) ZtJ3,1), (DELTACI,JsK), K=1-8) 
56 CONTINUE 
57 CONTINUE 
DO 59 1=l»4 
WRITE (6,260) THETA(I) 
WRITE (6950) (Y(1LyL), L=1y8) 
WRITE (6,55) 
DO 58 J=1,8 
WRITE(6,120) ZlJel)s (GAMMA I oJ eK), K=158) 
58 CONTINUE 
59 CONTINUE 
DO 62 I=1ls4 
WRITE(69270) THETA(T) 
WRITE (650) CYC LeL), L=1,8) 
WRITE (6,55) 
D0 61 J=1,48 


spl 


WRITE (62150) Z(Jol), (TCOMP (Ist eK), K=1,3) 
61 CONTINUE 
62 CONT INUE 
DO 64 I=1l+4 
WRITE(6,280) THETA(I) 
WRITE( 6250) (YC LoL), L=1,8) 
WRITE (6¢55) 
DO°63 y=1,8 
WRITE(69150) ZlJol)dsy (RCOMP(IoJ9K), K=1,8) 
63 CCNTINUE 
64 CONTINUE 
DO 66 I=1ls4 
WRITE(6e290) THETA(I) 
WRITE (6¢50) (YC 1leL)— L=1l,8) 
wRITE (6555) 
00 65 J=1,8 
WRITE (6,60) Z( deh)» (PTX(19JSeK)y K=1 98) 
65 CONTINUE 
66 CONT INUE 
DO 68 Il=ls4 
WRITE (6,310) THETA(TI) 
WRITE (6950) (Y¥(1leL)s L=1,8) 
WRITE (6,55) 
DO 67 J=1,48 
WRITE (6,60) ZO Jol), (PTTCI1eJSeK) » K=1,3) 


a2 


67 
68 


(ie 
10 
20 
30 
40 


50 


os) 


60 


TO 


80 


90 
L0G 


110 


CONTINUE 
CONTINUE 
WRITE (6,330) 
WRITE (6550) (Y€19L), L=198) 
WRITE (6,55) 
BO 72 J=ls8 
WRITE(6,60) ZlJel)»y (RIG SeK), K=1y8) 
CONTINUE 
FORMAT( F6.3) 
FORMAT( 2F6. 3) 
FORMAT(4F10.3) 
FORMATC£1",*UCTyJ9K) FOR SWIRL VANE BLADE ANGLE'»F5e1,! 
to'/) 
FORMAT(9Xs'Y' »8FS8 03) 
FORMAT(£O%,3X,'Z!) 
FORMAT ( F603 94X » BF Be2) 
FORMAT(*1*s*W(IyJeK) FOR SWIRL VANE BLADE ANGLE’ »F5.1,! 
+.'/) 
FORMAT("1",*RY( I, JeK) FOR SWIRL VANE BLADE ANGLE',F5.1, 
tots) 
FORMAT ( F6.394X» BF 83) 
FORMATC 1", *RZ(1,d,K) FIR SWIRLVANE BLADE ANGLE'»F5.1,! 
+.'/) 
FORMAT(*1*, *ALFAC I,Jd,K) FOR SWIRL VANE BLADE ANGLE',F5. 
+lyt.'/) 


Spe 


120 FORMAT(F6.354X,8F 8.0) 

130 FORMAT(®1%,*BETA( IyJeK) FOR SWIRL VANE BLADE ANGLE'*,FS5. 
$le*2 7) 

140 FORMATO 1% s*CCIsJdJeK) FOR SWIRL VANE BLADE ANGLE*',F5.1,° 
+.'/) 

150 FORMAT (F639 4X, 8F 8.3) 

160 FORMAT(.4 1%, *XCOMP(TyJyK) FOR SWIRL VANE BLADE ANGLE*,F5 


tely*.'/) 

170 FORMAT( *®1%,*RTOT(CI,JdJeK) FOR SWIRL VANE BLADE ANGLE',F5. 
+1,'.'/) 

180 FORMAT(*1*,*PHI([eJsK) FOR SWIRL VANE BLADE ANGLE*,F5.1 
tyf.t/) 

190 FORMAT(C*® 1*,*Al(IyeJeK) FOR SWIRL VANE BLADE ANGLE'*,F5.1, 
#*,9/) 

200 FORMAT( £1", *A2(IyJsK) FOR SWIRL VANE BLADE ANGLE'*,F5.1, 
+*.*7/)) 

210 FORMAT(*1%,*B(IeJsK) FOR SWIRL VANE BLADE ANGLE*,F5.1,! 
ata s, 

220 FORMATO * Ll 8 es *RYZ(CIeJen) FORSWIRL VANE BLADE ANGLE*,F5e.l1, 
+, 7) 

230 FORMAT(*1*s*RXZ0Ee9tJeK) FOR SWIRL VANE BLADE ANGLE® ,F5.1 
tes Sf) 

240 FORMAT(C 8 l* s*RXYC IT eJesK) FOR SWIRL VANE BLADE ANGLE',F5.1 
tees) 


250 FCORMAT('L§ s SDELTAC I,» 5eK) FOR SWIRL VANE BLADE ANGLE*,F5 
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eis 3 */) 

260 FORMAT('1®,*GAMMA(IyJ,/K) FOR SWIRL VANE BLADE ANGLE',F5 
els * .*/) 

270 FORMAT(*1*,*TCOMP(I,J+K) FORSWIRL VANE BLADE ANGLE',F5. 
mie. 74) 

280 FORMAT(*'1L*,*RCOMP(I,J+K) FORSWIRL VANE BLADE ANGLE',ES. 
tlyt.t/) 

290 FORMAT(*1*,*PTX(IyJ9K) FOR SWIRL VANE 3LADE ANGLE',F5.1 
te"et/) 

310 FORMAT(*L*,*PTT(LydeK) FOR SWIRL VANE BLADE ANGLE',F5.1 
tetet/) 

STIP 
330 FORMATC*L",*RL(J,K) FOR ALL EXPERIMENTS*/) 
400  FORMAT(8F8.4) 


410 FORMAT(*1*",* UMAX(1T)2 %54(F6ec9 65% 92X)) 
411 FORMAT(//,° WMAXCT)2 £54(F6.25'4'52X)) 
412 FORMAT(//,° AWUT): %,4(F6.29%9'52X)) 
413 FORMAT(//,° AU(T)2 %,40F6e29"%5' 52X)) 
414 FORMAT(//,° SAVOT)S ©54(F6e35%5 %e2X)) 
415 FORMAT(//;5°* SMAX(T) 2 %540F6.35%5'92X)) 
END 
SENTRY 


UMAX(IT): 475.00, 616-81, 759239, 820237; 


252 


AW( I): 


AU(T): 


SAV(T) > 


SMAX( 1): 


0.00, 


434254 


0.000, 


0.000, 


243-41, 


509. Ol,» 


Oe 478, 


02335, 


345.32, 


560-96, 


0.616, 


0.649, 


U( IlyJ9K) FOR SWIRL VANE 8LADE 


0.000 
0.050 
0.100 
0-150 
Je 200 
0-250 
0.300 
0-350 


Y 


0.00 


0-00 
0.00 
0.00 
50329 
554.37 
587.41 
540.67 
424%. 3% 


0.05 


0.00 
278.30 
376266 
524-27 
572.62 
515.44 
571-90 
416.06 


W(IeJeK) FOR SWIRL VANE 


0.000 
0.050 


O20 
0.00 


0.00 
142.03 


0.10 


398.18 
460-38 
492-4 84% 
532.85 
564.35 
a9 1. 92 
570.99 
416.06 


ZI0e tts 


462.45, 


0.640, 


0.544, 


ANGLE 15.0. 


585.08 
BT Seon, 
Sloot 
521.46 
587.42 
582.51 
527.95 
378.04 


595 249 
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